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EVOLUTION  OF  A  LOCALIZED  THERMAL  EXPLOSION  IN  A  REACTIVE  BA8 


T.  L.  Jackson*,  A.  K.  Kapila*  and  D.  S.  Stewart** 


Experimental  observations  of  ignition  in  premixed  gaseous 
reactants  indicate  that  perfectly  homogeneous  initiation  is 
practically  unrealizable.  Instead,  combustion  first  sets  in,  as  a 
rule,  at  small,  discrete  sites  where  inherent  inhomogeneities  cause 
chemical  activity  to  proceed  pref erential ly  and  lead  to  localized 
explosions.  Combustion  waves  propagating  away  from  these  "hot 
spots"  or  "reaction  centers"  eventually  envelope  the  remaining 
bul  k . 

This  study  examines  the  spatial  structure  and  temporal  evolution 
of  a  hot  spot  for  a  model  involving  Arrhenius  kinetics.  The  hot 
spot,  characterized  by  peaks  in  pressure  and  temperature  with 
little  diminution  in  local  density,  is  shown  to  have  one  of  two 
possible  self -similar  structures.  The  analysis  employs  a 
combination  of  asymptotics  and  numerics,  and  terminates  when 
pressure  and  temperature  in  the  explosion  have  peaked. 
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1.  Introduction 


This  paper  descibes  mathematically  the  birth  and  growth  of  a  hot 
spot,  or  localized  thermal  explosion,  in  a  premixed  reactive  gas. 

Experimental  observations  on  the  initiation  of  combustion  in 
gases  at  high  temperatures,  in  shock  tubes  and  elsewhere,  have 
demonstrated  conclusively  that  spatially  homogeneous  combustion  is 
essentially  an  unattainable  ideal.  In  fact,  ignition  first  sets  in 
locally,  in  small  volume  elements  at  discrete  sites,  where  chemical 
reaction  proceeds  pref erent i al 1 y  due  to  inherent  imperfections  in 
the  system.  In  due  course,  combustion  waves  originating  from 
localized  explosions  occuring  at  these  "hot  spots",  "reaction 
centers"  or  "exothermic  centers"  envelope  the  entire  reacting  mass. 
The  role  played  by  these  sites  as  precursors  of  more  dramatic 
combustion  phenomena  is  revealed  with  unsurpassed  clarity  in  Urtiew 
and  Oppenheim's  Cl  3  photographic  records  of  def lagration-to- 
detonation  transition  in  a  Hydrogen— Oxygen  mixture  confined  to  a 
tube.  These  photographs  show  that  as  the  deflagration  travels  down 
the  tube,  it  accelerates  and  evolves  into  a  highly  folded  turbulent 
flame,  preceded  by  a  so-called  precursor  shock.  Eventually,  an 
exothermic  center  is  formed  in  the  vicinity  of  the  flame,  near  the 
tube  wall.  The  localized  explosion  in  this  center  creates  a  blast 
wave  which  propagates  through  the  preconditioned  mixture  behind  the 
precursor  shock  and  ultimately  evolves  into  a  fully-developed 
detonation.  The  same  feature  appears  in  other  modes  of  detonation- 
initiation,  as  well  as  in  other  geometric  configurations. 

The  early  analyses  of  reaction-center  dynamics  are  due  to  Zajac 
and  Oppenheim  C23  and  Meyer  and  Oppenheim  C31.  In  these  studies 


the  reaction  center  is  assumed  to  be  a  spatially  homogeneous  source 
of  chemical  energy,  capable  of  expansion  and  separated  from  its 
inert  surroundings  by  an  impermeable  barrier,  across  which  only 
momentum  transfer  can  occur.  Either  by  prescribing  a  specific 
reaction  scheme,  or  by  specifying  an  energy  release  profile  within 
the  center,  the  above  authors  were  able  to  compute  the  resulting 
pressure  pulse. 

In  this  paper  the  reaction  center  is  treated  as  part  and 
parcel  of  the  reacting  medium  rather  than  an  isolated  entity  in  an 
inert  atmosphere,  and  is  found  to  have  a  definite  spatial 
structure.  The  aim  of  this  paper  is  to  describe  this  structure  and 
to  study  its  temporal  evolution  in  a  plane,  one-dimensional 
framework,  under  the  asssumption  that  the  reactive  gas  undergoes  a 
single,  one-step,  first-order,  irreversible  chemical  reaction  of 
the  Arrhenius  type.  One  may  argue  that  the  simple  overall  kinetic 
scheme  adopted  here  is  too  idealized  to  be  realistic.  However,  for 
large  activation  energies,  the  kinetics  does  capture  an  essential 
attribute  of  most  combustion  systems,  namely,  a  reaction  rate  which 
accelerates  rapidly  with  increase  in  temperature.  Thus  the  model 
is  quite  appropriate  for  studying  problems,  such  as  the  one  at 
hand,  which  owe  their  genesis  to  the  interaction  between 
gasdynamics  and  chemical  heat  release  at  highly  temperature- 
sensitive  rates. 

The  configuration  of  the  system  is  so  prescribed  as  to  provoke 
the  development  of  a  single  hot  spot,  and  this  can  be  accomplished 
in  a  variety  of  ways.  For  example,  the  shock-induced  thermal - 
runaway  studies  of  Clarke  and  Cant  E4D  and  Jackson  and  Kapila  C53 


Wfmwwmwrowwr 


-3- 

considered  a  semi-infinite  expanse  of  gas  ignited  by  a  piston- 
driven  shock,  thereby  creating  a  hot  spot  at  the  piston  face. 
Instead,  the  present  work  assumes  that  the  gas  is  confined  between 
two  parallel  planes,  and  that  its  initial  state  possesses  a  slight 
spatial  nonuniformity.  (In  a  practical  situation  these  non¬ 
uniformities  may  be  caused  by  a  variety  of  factors,  such  as 
turbulence,  interacting  pressure  waves,  or,  in  the  case  of 
condensed  explosives,  material  imperfections. )  The  mathematical 
model  leads  to  an  initial -boundary  value  problem  for  the  equations 
of  reactive  gasdynamics.  An  asymptotic  solution  is  developed  in 
the  limit  of  large  activation  energy,  and  the  analysis  is  carried 
as  far  as  the  end  of  the  localized  explosion  within  the  center. 

The  subsequent  expansion  of  the  center,  and  the  eventual  generation 
of  a  blast  wave,  will  be  the  subject  of  a  future  publication. 

The  temporal  evolution  of  the  explosion  occurs  in  two  stages, 
beginning  with  the  induction  stage.  Here  the  state  of  the  gas  is  a 
small  perturbation  of  the  initial  state  and  the  underlying  physical 
processes  are  those  of  linearized  acoustics  coupled  to  a  weak  but 
nonlinear  chemical  reaction.  The  reduced  equations  require  a 
numerical  solution  (see  C41  and  C53)  which  exhibits  local  thermal 
runaway.  Induction  is  followed  by  the  explosion  stage,  which 
consists  of  several  distinct  spatial  zones.  There  is  the 
practically  frozen  outer  zone,  and  a  rapidly  shrinking  inner  zone 
or  layer  in  which  intense  chemical  activity  leads  to  an  explosive 
growth  of  temperature  and  pressure.  Nonlinear  chemistry  is  again 
coupled  to  linearized  gasdynamics,  but  now  the  linearization  is 
about  an  atmosphere  undergoing  a  spatially  homogeneous  thermal 
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explosion.  As  the  layer  shrinks,  it  recedes  away  -from  the  outer 
zone,  thereby  creating  an  intermediate  zone  which  is  frozen  in 
time.  Although  highly  nonlinear,  the  explosion  stage  is  amenable 
to  analysis  because  gasdynamics  is  of  secondary  importance; 
temporal  variations  are  much  too  rapid  for  the  gas  to  undergo 
significant  expansion. 

For  the  specific  reaction  scheme  under  consi derati on  it  is  found 
that  the  reaction  center  can  have  one  of  two  posible  spatial 
structures,  depending  upon  whether  the  temperature  profile  within 
the  hot  spot  has  a  sharp  peak  or  a  rounded  peak  (Figure  1).  The 
former  typifies  hot  spots  originating  at  boundaries  (e.g.,  a  piston 
face),  and  the  latter  those  occuring  in  the  interior  of  the  vessel. 
These  structures,  which  will  be  referred  to  as  the  "Type  BH 
(boundary-type)  or  "Type  I"  (internal  type),  are  both  self-simi lar . 
The  former  is  described  below  in  detail,  with  only  the  results  for 
the  latter  given  in  section  6.  In  addition  to  these  two  structures 
there  exists  a  third,  described  briefly  in  the  Appendix;  it  is 
singular  and  corresponds  to  very  special  initial  conditions. 

The  specific  conf i gurat i on  under  study  here  was  also  examined, 
with  similar  methods,  by  Poland  and  Kassoy  C63.  Their  analysis 
differs  from  ours  in  one  crucial  respect;  they  considered  the 
distinguished  limit  in  which  the  spatially  homogeneous  induction 
time  at  the  initial  state  and  the  conduction  time  across  the  vessel 
are  of  the  same  order ,  i.e.,  the  Frank-Kamanetski i  number  S  i  s  of 
order  unity,  albeit  supercr i ti cal .  In  our  analysis  the  induction 
time  is  comparable  to  the  acoustic  time  across  the  vessel,  i.e.,  fi 
is  very  large.  In  physical  terms,  the  explosive  mixture  being 


considered  here  has  a  faster  reaction  rate. 


2.  The  Basic  Equations,  and  Bat up 

The  equations  of  reactive  gasdynamics  for  plane,  one-dimensional, 
unsteady  motion  are  C73 
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Here  p,  p,  T,  u  and  Y  are,  respectively,  the  gas  pressure,  density, 
temperature,  velocity  and  reactant  mass  fraction.  The  variables 
have  been  made  dimensionless  with  respect  to  a  constant  reference 
state  Pq,  Pq,  Tq  and  Yq.  Velocity  is  referred  to  the  acoustic 
speed  Cq,  defined  by 


co  =  LyPo/poll/2' 


time  to  tQ,  the  homogeneous  induction  time  at  the  reference  state, 
and  length  to  CQt  q.  The  diffusion  terms  have  been  left  out  because 
they  are  much  too  small  to  play  a  role  in  the  problem  under  study. 
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The  dimensionless  parameters  appearing  above  are  the  specific-heats 
ratio  y,  the  chemical  heat  release  ft  and  the  activation  energy  9. 

Let  the  reacting  gas  be  confined  to  the  interval  0  <  x  <  a.  At 
the  walls  the  appropriate  boundary  conditions  are 

(2.2)  u(0,t)  =  u (a , t )  =  0. 

The  initial  state  of  the  gas  is  taken  to  be  an  O(0-*)  perturbation 
of  the  spatially  homogeneous  and  stationary  reference  state,  i.e., 

(2.3a)  u (x ,0)  =  ©~1u1(x,0), 

(2.3b)  $(x,0)  =  1  +  0-1f1(x,O>  for  *  =  T,  p,  Y  and  p, 

where  the  precise  specification  of  Uj(x,0)  and  #^(x,0)  must  await 
the  next  section.  Note  that 

(2.3c)  Py (x , O)  =  pj(x,0)  -  Tj(x,0) 

in  accordance  with  the  gas  law  (2.1e).  An  asymptotic  solution  of  the 
i ni ti al -boundary-val ue  problem  (2.1)-(2.3)  is  sought  in  the  limit 
0  ■*  m,  with  ft  and  Y  fixed  and  0(1),  until  the  localized  explosion 
has  reached  completion.  The  various  stages  of  evolution  are 
detailed  in  the  following  sections. 

3.  The  Induction  Stage 

The  initial  conditions  (2.3)  suggest  that,  at  least  initially, 
the  state  of  the  gas  remains  an  0(©‘_1)  perturbation  of  the 


reference  state.  During  this  period,  referred  to  as  the  induction 
stage,  one  therefore  seeks  the  expansions 

(3.1)  u  -V  6  ^Uj  +  ■  .  •  ,  J'-l  +  0^§j  +  ..., 

for  #=  T,  p,  Y  and  P 

which,  upon  substitution  into  the  set  (2.1)  yield  the  leadi ng-order 
disturbance  equations 

(3.2a)  ( 3/ 3t  ±  3/3x)(p1  ±  /Uj)  *  /expiTj), 

(3.2b)  S/StC/Tj  -  (y-DpjD  ■  /expdj), 

(3.2c)  Pj  =  pj  -  Tj ,  3Yj/3t  =  -(l//9)exp(Tj)  . 

Except  for  the  nonlinear  source  term,  eqns.  (3.2a,b>  are  simply 
those  of  linearized  acoustics  in  a  uniform  atmosphere.  It  is  a 
simple  matter  to  integrate  them  along  the  characteri sti cs ,  as  was 
done  in  C43  and  C51  for  a  different  configuration.  During 
induction  it  is  enough  to  concentrate  on  the  variables  Tj,  pj  and 
Uj,  because  once  they  are  known,  the  first  eqn.  in  (3.2c>  yields  p ^ 
while  the  second,  combined  with  (3.2b)  and  integrated,  determines 
Yj  according  to  the  expression 

(3. 2d)  Tj  -  c  (y-n/yipj  +  =  t1(x,o)  -  c  (y-n/yipj  (x,o> 

+  /3Yj  (x  ,0)  . 

Equations  (3.2)  need  to  be  solved  numerically,  and  this  was  done 


and  -for  a  variety  of  smooth  initial  conditions  and  interval  lengths 
a.  A  high-resolution,  adaptive  DDE  integrator  was  employed  to 
integrate  along  the  characteristics.  All  computations  displayed 
thermal  runaway,  characterized  by  the  unboundedness  of  Tj  and  pj 
somewhere  in  the  interval  CO,al  at  a  finite  time  te-  The  numerical 
results  can  all  be  summarized  by  considering  two  representati ve 
cases,  for  which  the  initial  values  of  pressure  and  mass  fraction 
correspond  to  those  at  the  reference  state  and  the  initial  velocity 
i  s  zero,  i  .  e.  , 

(3.3a)  Pj(x,0)  =  Yj (x,0)  =  Uj (x ,0)  =  O, 

while  the  initial  temperature  perturbati ons  are  prescribed  as 

(3.3b)  Tj(x,0)  =■  bCl—  (x/a)  1  for  case  I,  bCl—  (x/a)^D  for  case  II. 

(Numerical  results  to  be  presented  below  correspond  to  a  “  0.B,  b  = 
0.5.)  In  both  cases  the  initial  disturbance  has  a  single  maximum 
at  x  *  0,  causing  it  to  become  the  site  of  thermal  runaway.  The 
essential  difference  between  the  two  cases  is  that  in  I  the 
temperature  disturbance  has  a  nonzero  spatial  gradient  (sharp  peak) 
and  in  II  a  zero  spatial  gradient  (rounded  peak),  at  x  =  0.  Thus  I 
typifies  a  hot  spot  located  at  the  boundary  (e.g.,  the  shock 
configuration  discussed  in  C41  and  C51),  and  II  an  internal  hot 
spot  (easily  visualized  by  a  symmetric  reflection  about  the 
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origin).  Henceforth  the  two  cases  will  be  referred  to, 
respectively,  as  Type  B  (boundary)  and  Type  I  (internal).  Their 
spatial  structures,  it  turns  out,  are  different. 

In  the  following  sections  the  Type-B  problem  is  discussed  in 
detail.  The  Type-I  problem  can  be  treated  analogously  and  is,  in 
fact,  slightly  simpler  to  analyze;  it  was  deemed  sufficient, 
therefore,  to  simply  state  its  solution  in  section  6. 

We  start  with  Figure  2,  which  displays  the  numerical  results  for 
the  Type-B  induction  solution.  The  four  graphs  there  exhibit, 
respectively,  the  profiles  of  Tj,  pj,  Uj  and  P^  against  x  for 
increasing  values  of  t,  upto  the  time  beyond  which  the  integration 
routine  was  unsuccessful  for  a  time  step  lO-^*,  thus  signalling  the 
imminence  of  blowup.  An  examination  of  the  Tj—  profile  near  blowup 
reveals  the  birth  of  a  boundary  layer  at  x  =  O.  Additional 
information  is  provided  by  Figure  3,  where  the  function  expC- 
Tj <0,t) 1  is  graphed  near  blowup.  The  straight-line  graph  in  the 
figure  has  slope  1.4  («/)  ,  and  a  t-intercept  equal  to  the  blowup 
time  te,  allowing  one  to  conclude  that 

(3.4)  Tj  (0,  t )  ~  —  JnC  /  ( te~ t )  3  +  o(l)  as  t  ->  te. 

Figure  4  displays  time  plots  of  the  solution  at  x  =0,  and  shows 
clearly  that  while  Tj(0,t)  and  pj(0,t)  become  unbounded,  p^Ojt) 
does  not.  Therefore,  pj(0,t)  must  have  precisely  the  same  leading- 
order  behavior  as  T^(0,t),  i.e., 

(3.5)  Pj^t)  ~  -£nCte-t]  +  0(1)  as  t  -*•  te. 
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To  summarize,  the  induction  stage  exhibits  the  classic  logarithmic 
singularity  at  spatially  homogeneous  thermal  runaway  E8D. 


4.  Type-B  Blowup  Structure 

Although  numerics  has  elucidated  the  temporal  character  o-f  the 
blowup  singularity,  -further  analysis  is  needed  to  ascertain  its 
spatial  structure.  This  will  be  done  by  examining  separately  the 
boundary  layer,  whose  emergence  has  already  been  noted,  and  the 
region  outside.  First,  it  is  convenient  to  introduce  a  new  time 
variable  r  via  the  expression 


(4.  1) 


r  =  t  -t,  r  >  0. 


Then,  following  simple  manipulations,  ears.  (3.2a,b)  transform  into 

(4.2a)  3pj/3x  -  TSuj/Sr  =  0, 

(4. 2b'  yauj/Sx  -  3pt/3r  =  (y-Dapj/Sr  -  yaTj/ar  =  yexptTj), 


where  the  dependent  variables  are  now  treated  as  functions  of  x  and 
r.  The  relevant  boundary  condition  is  the  first  of  (2.2),  re¬ 
written  as 


(4.3)  uj(0,r)  =  0. 

Elementary  manipulations  on  (4.2a,b)  and  (4.3)  yield  the  following 
integral,  which  will  prove  to  be  of  value  later  on: 
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te 

=  3q  exp  exp(Tj(0,t))dt], 

T 

where  3(r)  denotes  the  disturbance  temperature  gradient  at  x  =  0, 
and  3q  its  initial  value,  i.e. , 

(4.5)  3(c)  =  [3Tj  (x,  r)/3x]x=0,  3Q  =  3(t#-> . 

Recall,  from  (3.3b),  that  3q  vanishes  for  type  I  but  is  negative 
for  type  B.  Then  (4.4)  shows  that  3(c)  =  0  (rounded  peak)  for  the 
former  and  decreases  monotonocal 1 y  to  -e  (sharp  peak  approaching  a 
cusp)  for  the  latter  as  r  -*■  0+. 

4. 1  The  Boundary  Layer 

Turning  now  to  the  asymptotic  analysis  near  blowup,  eqns.  (4.2a,b) 
govern  the  region  outside  the  boundary  layer,  where  the  outer  limit 
process 


x  >  0  and  fixed,  c  -*■  0 

applies.  The  boundary  layer,  on  the  other  hand,  corresponds  to  the 
inner  limit  process 

s  >  0  and  fixed,  r  -*•  0, 

where  s(x,c)  is  the  spatial  coordinate  in  the  boundary  layer, 
reflecting  its  self-similar  structure.  The  shrinking  nature  of  the 
layer  requires  x  to  vanish  under  the  inner  limiting  process,  and 
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then  a  moment's  reflection  suggests  the  definition 


(4.6) 


x  /  r , 


which  assigns  coequal  importance  to  the  x-  and  r-deri vati ves , 
thereby  providing  the  richest  equations  for  the  inner  limit.  (It 
will  transpire  that  this  scaling  does  not  quite  cover  the  entire 
boundary  layer,  but  more  about  that  later.)  These  equations, 
obtained  by  transforming  (4.2a,b)  to  the  (s,r)  variables,  are 


(4.7a) 

(s3/3s  - 

r  3/  3r)  C  XT  j 

(4.7b) 

(s3/3s  - 

r 3/ 3r) p  j  + 

(4.7c) 

(s3/3s  - 

r3/3r)uj  + 

It  is  convenient  to  isolate  the  temporal  singularity  from  the 
spatial  structure,  by  setting 


(4.8a) 

T1 

=  —in ( /r) 

(4.8b) 

Pi 

=  -in(B1r) 

(4.8c) 

U1 

=  h (s, r) , 

where  the  yet  unknown  constant  Bj  represents  a  weak  influence  of 
the  initial  conditions  on  the  self-similar  boundary  1 ayer ,  and  will 
be  determined  in  due  course  by  matching.  The  structure  functions 
f,  g  and  h  are  assumed  to  be  o(l)  in  the  limit  r  ■*  0.  Substitution 
of  (4.8)  into  (4.7)  yields  the  structure  equations 
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(4. 

,  9a  > 

sfs  ' 

+ 

(y-i)hs 

=  e 

(4. 

,9b) 

s(fs 

-  gs>  - 

r(fr  - 

9r> 

(4. 

,9c) 

shs  - 

rhr  + 

(l/y)gs 

=  0 

The  only  boundary  condition  appropriate  -for  the  above  set  is  the 

wall  condition 

(4.  10)  h  (0,  r )  =  0. 

In  addition,  since  the  initial  data  are  smooth,  the  structure 
■functions  and  their  s-deri vati ves  are  required  to  be  regular  in 
Consider  the  asymptotic  expansions 


(4.11) 


$  ~  OjtrOijts)  +  a?(r)|?(sl  ...  for  $  =  f,  g,  and  h. 


^2 ' 1  •  *2 


as  r  O.  The  gauge  sequence  fon<r)>  is  not  yet  specified,  but  a 
clue  as  to  its  identity  is  provided  by  the  integral  relation  (4.4), 
rewritten  as 

T 

(4.12)  f  s  (0,  r)  =  30  r  exp  [-J*  £1/ (  ytlexpCf  (0,t>  >dt] 

*e 

in  view  of  the  scaling  (4.6)  and  the  prescription  (4.8).  For  small 
r  the  f-expansion  in  (4.11)  allows  the  above  relation  to  be  reduced 
further  to  the  asymptotic  form 


(4. 13) 


CTjfj'(O)  +  ...  = 


30rXexp[0(l)  -  £f1(0)//>J  t_1CT1(t)  dt  +...], 

0 


where 


Recall,  from  (4.5),  that  the  constant  3q  is  nonzero  for  the  Type-B 
problem.  Then,  the  assumption  that  fj'(O)  is  nonvanishing 
(involving  no  loss  of  generality)  leads  to  the  conclusion 

(4.15)  <r1(r)  =  rX 

if  the  two  sides  of  (4.13)  are  to  balance  at  leading  order.  With 
CTj  determined,  it  can  be  shown  that  the  expansions  (4.11)  proceed 
in  powers  of  rX. 

The  boundary-1 ayer  analysis  can  now  be  carried  out,  and  as 
hinted  earlier,  the  layer  is  found  to  have  a  two-sublayer 
structure.  It  is  convenient  to  refer  to  Figure  5  in  which  the 
various  spatial  regimes  near  and  beyond  blowup  are  displayed 
schematically.  OR  refers  to  the  outer  region  and  BL  to  the 
boundary  layer;  the  latter  is  subdivided  further  into  an  interior 
sublayer  Lj  and  an  exterior  sublayer  L^.  We  shall  first  examine 
the  interior  sublayer,  show  that  it  becomes  nonuniform  for  large  s, 
determine  the  appropriate  scaling  and  expansions  for  the  exterior 
sublayer,  and  demonstrate  that  the  latter  merges  smoothly  into  the 
outer  region.  Only  one  or  two  terms  of  the  expansions  in  each 
region  will  be  computed;  continuation  to  higher  orders  is 
straightforward  though  increasingly  complex  al gebrai cal  1 y . 

4.1.1  The  Interior  Sublayer  Lj 

Substitution  of  (4.11)  into  (4.9)  yields  the  leading-order 


structure  equations  for  the  inner  sublayer. 


(4. 16a) 

ill 

■4* 

-  (X+Dfj  +  (M)hj'  = 

(4. 16b) 

s(f  l 

-gi')  -  -  hj' 

(4. 16c) 

in 

T 

h* 

-  Xh j  +  (1  /y)gj  '  =  0. 

The  boundary  condition 

(4. 16d )  h j (0)  =  0 

comes  -from  (4.10),  and  the  solution  is  restricted  additionally  by  the 
requirement  that  it  be  regular.  If  g^  and  hj  are  eliminated  -from 
(4.16a~c),  the  result  is  the  third-order  equation 

(4.17)  C (s3-s) f  j '  +  AC (l-3s2)f  j ' 3 '  -  (s2-l / /) f t ' ' 

+  (A+l) (3A-4)sf j '  +  (X2-l) (2-A)f j  -  O 


for  fj.  The  points  s=0  and  s=l  are  singular  points  of  this 
equation  and  the  three  linearly  independent  solutions  have  the 
asymptotic  behavior 


1,  s  and  sins  as  s  -»  0,  and 

1,  1-s  and  1 1-s  I  ^  as  s  -*•  1. 


In  general  one  can  expect  a  one-parameter  family  of  regular 
solutions  to  exist,  and  numerical  computations  verify  that  such  is 
indeed  the  case.  A  convenient  parameter  is 


(4. 18) 


Aj  =  f  j (O) • 


With  known,  h^'  can  be  eliminated  from  (4.16a,b)  to  obtain  a 


-first-order  di  f  f  erenti  al  equation  for  g.  whose  regular  solution 


turns  out  to  be 


(4. 19a) 


gt  *  C  <y+l/A)/(y-l)  If  J  -  CsX/X(y-l)  x  xf1(x)dx, 

0 


and  then,  (4.16c)  integrates  to  give 


(4.19b)  h1  =  g1'/(/X)  -  CsX/0'X)lJ‘  x"Xg1  '  '  (x  )  dx  , 

0 


where  regularity  has  been  imposed  again.  Thus  the  full  solution  at 


this  order  depends  on  the  single  parameter  Aj.  Graphs  of  f j ,  gj 


and  hj  for  Aj  *=  1  are  drawn  in  Figure  6. 


At  this  stage  the  solution  (4.8)  has  the  following  expansions 


in  the  interior  sublayer: 


(4. 20a)  Tj  ~  -£n(/r)  +  rAf1(s)  +  ..., 


(4.20b)  Pj  ~  -jen(Bjr)  +  rxgj (s)  +  ..., 
(4.20c)  ut  ~  rXh1(s)  +  ...  . 


In  order  to  determine  the  spatial  extent  of  Lj  one  needs  the 


asymptotic  behavior  of  fj,  g^  and  hj  for  large  s.  This  is  easily 


obtained  from  (4. 17)  and  (4. 19) ,  as 


(4.21a)  f1  ~  -A1aCs<2>'  l)//  +  s  1//(Af£ns  +  Bf ) 
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(4.21b) 


(4.21c) 


+  CrS 


-(y+n/y 


n  ~  +  D_s<y-i>/y  +  s-i/y^ns  +  B  > 


+  Cg.“(y+i)/y  + 


3, 


hj  -v  -AjaCs(y  ^^/y(Ahin«  +  B^)  +  C^s 


-i/y 


.  1. 


Here  <x,  and  Cf  are  constants  with  values 


a  =  2.660,  Bf  =  -0.236,  C*  =  0.0737 


obtained  by  integrating  the  -f  j  equation  (4.17)  numerically.  The 
remaining  constants  appearing  above  are  given  by 

Af  =  -(2y-i) (y-i)2/(2y3) ,  Ag  «  (y+n At/ (y-i> , 

Ah  -  2Afy/(y-i)2,  Bg  ■  Af/(y-i)  +  (y+i)B^/(y-i) , 

(4.22)  Bh  =  2Bfy/(y-i)2  -  y(3y-i>Af/(y-i)3, 

cg  =  (2y+i)Cf/(2y-2) ,  ch  =  -3ycf/(y-n, 

Dg  -  -3cf y3/ (y-i)2. 

The  range  o-f  validity  of  the  expansions  (4.20)  can  now  be  determined. 
For  example,  substitution  of  (4.21a)  into  (4.20a)  suggests  that  the 
latter  becomes  nonuniform  when 


rxB(2y-i)  /y 


0(1) ,  i .e. , 


0(r~w)  , 


where 


(4.23)  u  =  (y-l)/(2y-l) 


% 
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and  the  definition  (4.14)  of  X  has  been  invoked.  Correspondingly, 


x  =  0(ry/(2/-1))  =  o(l). 


The  smallness  of  x  indicates  that  although  one  has  reached  the  edge 
of  Lj,  the  outer  region  is  still  too  far.  The  need  for  an  exterior 
sublayer  is  therefore  apparent. 


4.1.2  The  Eh ter lor  Sublayer  Lg 

In  this  sublayer  the  appropriate  variables  are  ?  and  r,  with  f 
defined  by 


(4.24) 


5  —  s  s  1*. 


The  expressions  (4.B)  for  T1?  p^  and  Uj  hold  again,  provided  f,  g 
and  h  are  now  treated  as  functions  of  ?  and  r.  The  structure 
equations,  obtained  from  (4.9)  by  transforming  from  s  to  $,  are 


(4.25a)  <l-w)ff5  -  rfr  +  h?  =  e*-l, 

(4.25b)  ( l-/i)  ?  (f  ?-ge>  ~  r(fr-gr)  -  h?  =  O, 


(4.25c)  (l-w>?h?  -  rhr  +  (l//)rw  g?  =  O. 


Matching  requirements  imposed  by  Lj,  obtained  by  substituting 
(4.21)  into  (4.20)  and  then  employing  (4.24),  are 

(4.26a)  f  ~  -Ajo:  f<2/-l)/y  +  0(r2/i  fnr), 

(4.26b)  g  -v  -Atot  ?<2/-l)/y  +  0(rw), 


y.<\ 

/V. 
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*  a  *  I 
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(4.26c)  h  a.  -Aj<x  vu  f ( y  1)/yCAh<-*i  Inr  +  £n$)  +  Bhl 

+  0(r2w)  ,  as  f  +  0. 

Guided  by  these  the  L^-solution  is  sought  in  the  -form 

(4.27a)  f  -v  F0(?)  +  .  .  .  , 

(4.27b)  g  Gq(  f )  +  .  .  .  , 

(4.27c)  h  -v  ruC£nr  H0(?)  +  H1(f)3  +  ...  . 

Substitution  into  (4.25)  leads  to  the  differential  equations 

(1-W)#F0'  *  exp(F0)-l, 

Fo'  -  Go'  =  °» 

(l-fj)(H0'  -  vH0  =  0, 

<l-V>fH1'  -  wHj  =  H0  -  (  1  / /) Gq  ' , 

whose  solutions,  subject  to  the  matching  requirements  (4.26) 

(4.28a)  F0(f)  =*  Gq  (  f )  =  -  £nC  l+A^f  <2y“1  > /y3  , 

(4.28b)  Hq(?)  =  -A1«C(y-l )/y2yf(y~1)/r, 

(4.28c)  H  j  ( f )  =  Aj«  ?<y~1)/y[<(2y-l) /y2}i£n(  +  F0($>>  -  B 

Thus  the  Lg--solution  can  be  written  as 

(4.29a)  Tt  -v  -Int/r)  +  FQ  ( ?  >  + 

(4.29b)  px  -v  -£n(B1r)  +  FQ  ( f )  + 

(4.29c)  u,  ^  rw  tint  Hft(()  +H,(f)l  +  ...  . 


-20- 


Qne  must  consider  the  behavior  of  this  solution  -for  large  S  in 
order  to  assess  the  spatial  extent  of  the  sublayer  l_£.  This  can  be 
done,  for  example,  by  substituting  the  large-f  behavior  of  (4.28a) 
into  (4.29a).  The  result  is  the  expansion 

Tt  -v  -C(2y-l)/y>ln(r>'/ (2/-1)  ?)  -  ln(A1«)')  +  ..., 
as  $  ■*  e, 

which  clearly  becomes  disordered  when  g  =  0 ( (2/-1 ) > . 
Correspondingly,  x  =*  0(1),  indicating  that  the  edge  o-f  the  boundary 
layer  has  now  been  reached.  The  next  step  is  to  see  if  the 
boundary  layer  merges  smoothly  with  the  region  OR. 

4.2  The  Outer  Region  OR 

In  the  outer  region,  where  x  and  r  are  the  proper  variables,  the 
solution  can  be  expanded  as 

(4.30)  ~  #jQ(x,te)  +  r  +  ‘  •  *ar  =  Tj,  pj  and  Uj, 

where  the  leading  terms  are  the  numerically  obtained  limiting 
values  at  blowup  and  the  higher-order  terms  can  be  computed  from 
(4.2a,b)  under  the  outer  limit  process.  It  is  a  strai ghtf orward 
matter  to  establish  that  a  match  of  (4.30)  with  the  l_£-solution 
(4.29a-c)  requires  the  following  asymptotic  behavior  of  the  outer 
solution  at  blowup: 

(4.31a)  Tj  ~  -{  (2/-1 )  /  /Jinx  -  «n(A1«)')  +  ..., 


w 

Si 

V 

s. 
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(4.31b)  -  -C  (2/-1) //>*nx  -  £nif\l*B1)  ♦ 

(4. 31c  >  Uj Aia  x<y_1/y)  [€(2y-l)/y2>C<(y-l)/y>lnx 

+  £n ( A^ «) 3+B^ j  +  .  ..,  as  x  O. 

A  careful  examination  of  the  numerical  solution  does,  indeed, 

con-firm  this  behavior.  The  constants  Aj  and  B^,  the  only  ones  yet 

undetermined,  can  then  be  -found  by  comparing  the  above  expansions 

with  the  numerical  solution.  The  comparison  is  made  at  the  "edge" 

y  /  ( 2  y- 1  > 

of  the  boundary  layer,  i.e.,  -for  (x,r)  satisfying  r  <<  1,  r 
<<  x  <<  1.  It  should  be  emphasized  that  the  structure  of  the 
blowup  singularity  is  influenced  by  the  initial  conditions  only  via 
these  constants;  otherwise,  the  solution  has  a  universal,  self¬ 
similar  structure. 


4.3  Summary 

The  near-blowup  analysis  is  now  complete,  and  can  be  summarized. 
In  the  interior  sublayer  Lj  the  expansions  are 

(4.32a)  T  ~  1  +  0~1C-fn(yr)  +  rXfj(s>  + _ 3  +  ..., 

(4.32b)  p  -v  1  +  ©~*C-£n(Bjr)  +  (s)  +  ...3  +  ..., 

(4.37.C)  u  -v  0-1CrN-»^(s)  +  .  ..]  +  ..., 

where  f  j,  gt  and  h^  are  defined  by  (4.17)  and  (4.19).  In  the 
exterior  sublayer,  the  solution  is 


(4.33a) 

(4.33b) 


T  ~  1  +  e“1C-fn(yc)  +  F0(f>  +  ...3  +  . 
p  ~  1  +  C-£n  (Bj  r)  +  Fq(?)  +  ...3 
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(4.33c)  u  9  lZzv  fnr  Hn(?)  +  ru  H,(£)  +  .  .  .  3 


’O’  c  1 

where  Fq,  Hq  and  are  given  by  (4.28).  In  the  outer  region  the 
expansions  take  the  form 


(4.34a) 


T  -v  1  +  Q  1CT10(x,te)  +  0  ( r )  3  +  .. 


(4.34b) 


p  -v  1  +  0  *Cpj0(x,tfi)  +  0(r)3  +  ..., 


(4.34c) 


u  -v  0  1Cu^0(x,te)  +  0(r)3  +  ..., 


where  Tjq,  p10  anc*  ul0  ara  terminal  values  of  the  induction 

solution,  determined  numerically. 

The  remaining  variables  p  and  Y  can  be  computed,  upto  0(<?~1),  by 
appealing  to  the  first  equation  of  (3.2c)  and  (3. 2d).  The  results 
are 


( 4 . 35a  > 
(4. 35b) 


P  ~  1  +  9  1C«n(y/B1)  +  (gj-fj)  + 


3  + 


Y  -v  1  +  (B/9y)“1C«n(b1r)  +  zX  f  ( /-l )  gj-Xf  + 


3  + 


in  L 


I  ’ 


(4. 36a) 
(4.36b) 


P  -  1  +  G  1 C In ( //B ! )  +  ...3  + 


Y  -v  1  +  (OfiY)  1C«n(b1r)  -  FQ  ( ? )  +  ...3  +  ... 


in  ,  and 


(4.37a) 
(4. 37b) 


P  -v  1  +  e  1Cp1Q(x  ,tp)-T10(x  ,te>  +  0  (  r )  3  +  . 
y  -v  i  +  (Bfl)_1C(:(y-i)/y>p10(x,te)-T10(x,te) 
+  b(l-x/a)  +  0(r)3  +  ...  . 


a 


•Av, 


w 

•/ v 


*  «r-  m 


4$Cr 


.  -v 


.--Si 


«/  V  *.* 


r\ 

A 


*v  **,  . 
V.v.v 


V. .  , 


m 


r. 


-23- 


in  OR.  The  constant  bj  appearing  in  (4.35b>  and  <4. 36b)  is  given 
by 

( 4 . 38 )  bj  =  expC  yb  +  yin Y  -  ( y-1 ) InBj 1 . 

Observe  that  the  BL-solutions  (4.32),  (4.33),  (4.35)  and  (4.36) 

break  down  when  -inr  =  0(9),  signalling  the  end  o-f  the  induction 
stage,  and  the  onset  o-f  explosion.  In  contrast  the  OR-sol ut i ons , 
(4.34)  and  (4.37),  su-f-fer  no  disordering  and  in  -fact,  become 
increasingly  accurate  as  r  ■»  0. 


...  .u».h • 
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5.  The  Type-B  Explosion  Stags 

The  nonuniformity  just  encountered  decrees  that  -further 
evolution  in  the  boundary  layer  occur  on  the  new  time  scale  cr, 
defined  by 


(5.  1) 


For  a  =  0(1)  the  limit  G  ■»  e  corresponds  to  a  time  interval  of 
exponential  brevity;  its  role  in  the  evolution  of  thermal 
explosions  was  first  recognized  and  exploited  by  Kassoy  C81.  The 
two  sublayers  comprising  the  boundary  layer  must  again  be  examined 
in  turn.  In  fact,  we  shall  find  that  as  the  boundary  layers 
continue  to  shrink,  an  expanding  void,  or  an  intermediate  region 
(denoted  by  IR  in  Figure  5),  is  created  between  the  sublayer  l_E  and 
the  outer  region  OR;  this  region  begs  a  separate  treatment. 


5.1  The  Interior  Sublayer  Lj 

The  spatial  coordinate  in  this  region  remains  s,  now  written  as 


s  =  x  /  r  s  e^CT  x  , 


thereby  expressing  explicitly  the  continuous  shrinkage  of  the 
region.  In  the  (s,a)  variables  eqns.  (2.1)  transform  into 


(5.  2a) 


(5.2b) 


G  1  pa  +  sps  +  (pu)s  =  0, 

p[0-1  Ua  +  5Uj]  +  ( l//)ps  +  Puug  =  o. 
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(5.2c) 


i(y-i) /yjza  1  pa  +  sps  3 


pie'1  Ta  +  sT*l  - 

+  uCpts  -  c(y-i)/y>psD  =  w, 

<5. 2d)  pee-1  +  sYs:  +  upYs  =  -(1/0>W, 

(5. 2e)  p  =  PT, 

where 

(5.3)  W  =  0— 1  pYexp C £>(  1— <j— 1 /T )  3. 

The  boundary  condition  (2.2)  is  rewritten  as 

(5.4)  u(0,a)  =  0. 

At  -fixed  s  the  solution  must  match  with  the  induction  zone  as  a 
O.  To  obtain  the  necessary  conditions  one  applies  the  "explosion 
limit"  a  fixed,  &  ->  ®  to  the  Lj-solution  (4.32),  (4.35)  and  gets 


(5.5a) 

T 

1+cr  -  Iny  +  . 

.  .  + 

fiCf  x  (s)  +  .  .  . 

(5.5b) 

P 

l+<7  -  8~ *  #nB^  +• 

...  + 

5Cgj (s)  +  . . 

•3, 

(5.5c) 

u 

•V 

fiChj (s)  +  . . . 1 , 

(5.  5d ) 

P 

1  +  0-1  Mn ( y/Bj ) 

+  ... 

+  SCg1(s)-f1 

(s)  +  . 

( 5 . 5e ) 

Y 

l-a/(fly)  +  (9f}y)~ 

*  £nb  ^ 

+  ... 

+ 

sifly)  *c  (y-i>g1(s>-yf  j 

<s) +  . . . 1 ,  as 

a  ■*  0. 

Here , 

(5.6)  6  =  e~X  e-cr0X, 


where  X  was  defined  in  (4.14)  and  b^  in  (4.38).  These  conditions 
reveal  that  spatial  variations  in  the  explosion  stage  appear  only 
at  the  (exponentially  small)  0(fi)  level,  thereby  suggesting  that 
the  solution  is  spatially  uniform  to  all  algebraic  orders  in  6.  In 
other  words,  the  structure  of  the  interior  sublayer  consists  of  an 
extremely  weak  chemico-acoustic  field  superimposed  over  a  uniformly 
exploding  atmosphere.  Accordingly  one  seeks  expansions  of  the  form 

(5.7a)  u  ~  8  Uj(s,(r)  +  ..., 

(5.7b)  $  -v  $q(ct;£?)  +  8  $^(s,a)  +  ...,  for  $  *  T,  p,  p,  and  Y, 

with  the  understandi ng  that  the  $q  contain  all  terms  of  algebraic 
orders.  Substitution  into  (5.2)  finds  the  $q  satisfying  the 
standard  equations  of  constant-volume  thermal  explosion  C83  ,  i.e., 

(5.8a)  3Pq/ 3a  =  0,  Pq  =  PqTq, 

(5.8b)  (1/X) P03T0/3a  =  -flP03Y0/3a  =  WQ  2  pQYQexp C G ( 1— a— 1 /Tq) 3 . 

The  solution,  subject  to  the  matching  conditions  (5.5),  is 
(5.9a)  §0  ~  §qq  +  9~ 1  +  ...,  for  $  =  T,  p,  Y  and  P, 

where 

(5.9b)  T00  =  p00  -  (1-a)-1,  Y00  =  (l  +  y/3-T00>/<jrfl)  ,  P00  -  1, 

(5.9c)  T01  =*=  -(l-a)_2«n [/(l-a)2Y00] ,  pQ1  =  inOVBj), 

(5. 9d )  p01  =  T01  +  p01/(l-a>,  pyv0l  =  C£n(b1//)-T013. 


The  structure  functions  satisfy,  to  leading  order,  the  equations 


(5. 10a) 
(5. 10b) 
(5- 10c) 
(5. lOd ) 
(5.  lOe) 


(s3/3s-X)p1  +  pq3 uj/3s  =  0, 
P0(s3/3s-X)u1  +  (l//)3p1/3s  =  0, 
(s3/3s-X)Cp0T1-C(y-l)/)'>p1:  =  TjWq/Tq2, 
P0<s3/3s-X) Yt  =  -  TjWQ/fflTQ2), 

Pi  _  PqT1  ~  T0^1  =  °> 


where  Wq  was  defined  in  <5. 8b).  Replacement  of  Tq  and  Wq  by  their 
1 eadi ng-order  values  from  (5.9),  followed  by  the  use  of  the  trans¬ 
formations 


(5.11)  Uj  —  (1 -a)1/2  Uj,  s  =  s(l-a)  *^2, 


reduces  the  set  (5. lO)  to 


(5. 12a) 

(s3/3s-X)T1 

+  (y-nsuj/ss 

(5. 12b) 

(s3/3s-X)  (T 

j-pj)  -  3u^/3s 

(5. 12c) 

(sS/Ss-Xluj 

+  (l/yispj/ss 

(5. 12d> 

Pi  =  <P!-Ti 

>^00* 

(5. 12e) 

(s3/3s-X) Yj 

+  Cl 

Eqns.  (5. 12a-c)  are  identical  to  (4. 16a-c)  if  Tj,  pj, 
the  former  are  identified,  respect i vel y ,  with  fj,  g^  , 
the  latter.  Following  the  arguments  of  section  4.1.1 


A  A 

u j  and  s  in 
h^  and  s  in 
,  therefore. 


one  is  led  to  the  solution 
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(5. 13a) 


(5. 13b) 


(5. 13c) 


=  CAj  (<r)  /(\t  3f  j  (s)  , 

=  CAj  (<r) /Aj  <s) 

=  (1-<t)1/2u1  =  ( 1-cr)  1/2CA1  (<x> /Aj  ]hj  (s)  , 


where  Aj  is  the  constant  introduced  earlier  in  (4.18).  The 

A 

amplitude  function  A^(ct)  is  unknown  at  this  stage,  and  will  be 
determined  by  matching  with  the  exterior  sublayer.  So  far  we  only 
know  its  initial  value  as  a  result  of  matching  with  the  induction 
solution  (4.32),  i.e.. 


(5. 14) 


At (0)  =  At. 


It  is  now  a  simple  matter  to  solve  (5. 12d)  for  p^  ,  and  compute  Yj  by 
integrating  (5. 12e)  subject  to  the  regularity  requirement.  The 
resulting  expressions  are 


(5.15a)  Pt  =  (1-a) CAj (a) /A13Cg1 (s)-f j (s) 3, 


(5. 15b) 


Yj  =  CAj  (a)  /  (flAj  )  3C  C  (  /-l ) //>g1  (s)~f  j  (s)  1. 


Both  the  spatially  uniform  and  the  spat i al 1 y— varyi ng  components 
of  the  expansions  (5.7)  are  thus  determined  at  leading  orders, 
although  the  latter  involve  Aj (a)  which  is  still  to  be  found.  It 
is  worth  noting  that  the  spatial  structure  of  ths  solution  is 
essentially  the  same  as  it  was  at  induction-stage  blowupf  the 
scalings  (S. 11)  simply  reflect  the  temporal  evolution  of  the 
acouat i c  speed . 

As  in  section  4.1.1,  the  Lj-solution  breaks  down  for  large  s. 
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m  . 
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•  .  ■  ,  ’  * 


.V-V;." 


*  »"V 


*" ,  ■"  •f'm  *«*.  "'k  *(*.  *  •  ,  *' , 
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IWWW 


the  nonuniformity  now  occuring  (see  the  expansion  (4.21))  at 
OS 5^2/-!) //  _  Q(i)_  One  is  then  led  to  the  exterior  sublayer. 


9.2  The  Em tar i or  Sub 1  ay ar  LE 

Here  the  proper  variables  are  f  and  a  where  £  is  now  related  to  x 

A 

and  s  via  the  expressions 


(5. 16) 


$  =  x 


,0o y / <2/-i ) 


(  l-CT) 


-1/2 


a  e"*^. 


'>■ 


and  u  was  defined  in  (4.23).  In  the  new  variables  the  full  equations 
(2.1)  read 


(5.17a)  fiT1  pa  +  (1 -*j)£p?  +  a~0aulpu) ?  =  O, 

(5.17b)  pe0ou  CO'1  ua  +  (l-w)tu53  +  (l//)p?  +  puu?  =  0, 
(5.17c)  ptO'1  Ta  +  (l-ii)fT?3  -  t(y-l)/y>CO_1  pa  +  (l-w)fpf3 
+  e~0au  uCpTg  -  {(M)//)p{]  =  W, 

(5. 17d)  PC©-1  Ya  +  (l-y)fY?I  +  e~0au  upYf  =  -<l/fl)W, 

(5.  17e)  p  =  PT, 


where  W  retains  the  definition  (5.3).  The  solution  is  subject  to 
the  following  matching  conditions  imposed  by  Lj: 


(5. 18a) 

T 

JV  Too  +  0  1 

CToi  - 

P£ (2/-1 ) //3  +  .... 

(5. 18b) 

P 

~  T00  +  0  1 

CT01  + 

p01/(i-<t>  -  P5<2y-i)/y3  + 

(5. 18c) 

u 

•V  -e-^W  p? 

( r- 1 )  /  / 

[Ahpa  ♦ 

9  1  CAh£nC < l-o) 1/2f>  +  Bhl]  +  ..., 
p  -v  1  +  9~l  PQl  + - , 


(5. 18d ) 
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(5. 18e> 


~  Y00  +  (Ofly)  ^Jntbj//)  -  T01  + 


P$ 


(2y-i )  // 


] 


.  .  .  ,  as  f  -*  O. 


Here, 

(5.19a)  P  (a)  =  aAj  (a)  ( 1-a)  (2y~1  )  /  (2/)  , 

wi  th 

<5. 19b)  P(O)  =  aAj , 

from  (5.14).  In  obtaining  the  conditions  (5.18)  Me  have  employed 
the  expansions  (4.21)  and  the  solution  (5.13);  the  variables  with 
double  subscripts  are  the  spatially  homogeneous  functions  appearing 
in  (5.9).  It  turns  out  that  compliance  with  these  conditions  also 
ensures  temporal  matching  with  the  induction  stage.  The  L^- 
solution  is  now  sought  in  the  form 


(5.20a) 

T 

T00  + 

0~i  T1(?,ct) 

+  •  1 

( 5 . 20b ) 

P 

•V 

+ 

o 

o 

I- 

0-1  pj  <  5  ,  <r) 

+  •  -  •  t 

( 5 .  20c ) 

u 

e-0au 

Cuq  <{,o)  + 

Uj ( # ,CT)  +  ...], 

(5. 20d ) 

Y 

*\r 

Y00  + 

e~l  Yj ( g ,ct) 

(5.  20e) 

P 

■*N# 

1  +  0' 

»-* 

t> 

o 

+ 

• 

■ 

• 

• 

Substitution  into  (5.17)  shows  that  (5.17a)  is  satisfied 
identically  to  0(8_1).  At  0(1),  (5.17b)  reduces  to 


Km  « -m - 


Ji  vVV AW  *  - »  -  *  ^ 


V. 

vjyi 

v 

S* 


^.V.v; 


,*  ,*  V  ' 

... 


*\ 

•  V  V. 


v* 


T.''r.'‘rri 


il-v) 53u0/3?  -  wu0  =  O 


whose  solution  subject  to  the  matching  requirement  (5- 18c)  is 

(5.21)  uQ  =  <(y-l)/y2>aP§(>'-1)/>'. 

At  O(0_1),  (5. 17e)  yields 

(5.22)  pi  ■  Tj  +  p01/(l-a), 

while  (5.17c)  reduces  to 

( i-<r) — 2  +  (i-«)  f  aTj/at  -  c  (y-i) /y>  [d-a)“2 

♦  ( l-#i)  ?  3pi  /  3?]  =  Y00  expt ( 1-ct) 2Tj 3 , 

and,  in  view  of  (5.22),  simplifies  further  to 

(5.23)  (l-a)~2  +  (l-«) f3T 3f  =  yYQ0  expf ( l-a> 2Tj 3 . 

Its  solution,  consistent  with  the  matching  condition  (5.18a), 

(5.24)  ?!  =  T01  -  ( 1  -a)  ~2  in  [l  +  ( 1  -a)  2P f  <2  Y~ 1 )  '  y]  . 

■I*  *V 

With  Tj  known,  (5.22)  defines  pj.  In  order  to  determine  Yj 
consider  (5.17d)  at  0<£)~*)j  it  yields 

-(fly)  _1  (  1-ct)  “2  +  (l-p)?3Y,/3f 
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=  (  — l/fl)Yoc  ex  p  C  ( 1  —  <y)  2T  j  D  . 

When  linearly  combined  with  (5.23)  the  above  equation  leads  to 

dLTl  *  fl/Yjl/af  =  0. 


The  matching  condition  (5.  18e)  then  provides  the  -following  expression 
tor  Y j  : 


(5.25) 


fl/Yj  =  #n(b1/>')  -  Tj 


It  now  remains  to  determine  u^,  and  the  -function  P(<r>  (or, 
equivalently,  A^(ct)).  Both  are  obtainable  from  (5.17b)  which,  at 
0 ( 9~ 1 ) ,  reads 


(5.26) 


(l-p>  Sauj/35  -  =  -3uq/3ct  -  (l//)3p1/35. 


With  Uq  and  p^  known  (see  (5.21),  (5.22)  and  (5.24)),  the 

general  solution  of  the  above  equation  can  be  written  as 


(5. 27) 


Uj  =  {  (2y-l  ) //>  e  (  y  1>/y[K  -  <  (  y-1  > /y2>  (trP)  ’  j«nf 
+  (P/X)  JnC$  (2/_1>  /y/d+P(  l-cr)2f  <2y_1 }  /y>  ]]  , 


where  K(ct)  is  the  integration  "constant”.  As  5  -*■  O,  uj  has  the 
asymptotic  behavior 


(5.28) 


ut  -  c  (2y-i  > //>  f  (  ¥  l)  /  y  [i-ay-u  /  y2y  (op)  ■ 


V 


C <2y-l ) /y2>PDlnf  +  k]  +  ... 


>-1 


which  must  agree  with  the  0(0  )  term  in  (5.18c).  Matching  the  Inf 

terms  yields  the  differential  equation 


CctP  (ct>  1  '  =  P 


whose  solution,  subject  to  (5.19b),  is 


(5.29a)  P  =  «Aj , 


m 

'■.yM 

•ray 


or,  equivalently, 

(5.29b)  =  A1(l-a)~<2y_1,/<2/). 

With  P  determined,  matching  of  the  f-i ndependent  terms  in  (5.28) 
and  (5.1Bc)  yields  K: 

(5.30)  K  =  (ocAj/n  C  (1/2)  In(l-a)  -  O'2/ (2/-1 )  }Bh3  , 

where  the  constant  B^  was  defined  in  eqn.  (4.22).  The  L^- 
solution  at  the  explosion  stage  is  thus  complete. 

It  is  instructive  to  compare  the  solutions  in  the  two  sublayers. 
In  each  the  background  field  is  that  of  a  spatially  homogeneous 
thermal  explosion,  but  the  superimposed  spat i al 1 y-vary i ng  field  is 
quite  different,  both  in  amplitude  and  structure.  In  Lj  the 
spatial  component  is  exponentially  small  in  amplitude  but  has  a 
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chemi co-acousti  c  character;  all  disturbances  to  the  background 
homogeneous  field  are  of  the  same  size.  In  the  spatial 
variations  in  T,  Y  and  p  are  O(0-*) ,  while  those  in  u  and  P  are 
exponentially  small,  i.e. ,  the  evolution  is  essentially  due  to 
constant-volume  chemical  amplification  of  a  spati al 1 y— non- 
uniform  field,  with  gasdynamics  playing  a  very  minor  role. 

As  a  increases,  T  and  p  increase  and  Y  decreases,  in  both  the 
sublayers.  Eventually,  p  and  T  peak  when  Y qq,  the  leading  term 
in  Y,  vanishes.  This  happens  at  (see  (5.9)) 

(5.31a)  a  =  (IY  /  (  1+/9 Y)  , 

and  the  peak  values  are 

(5.31b)  T  -  1+/3/,  p  ~  1  +f)Y. 

At  the  same  time,  the  O(0-1)  term  (in  T,  say:  see  (5.20a),  (5.24) 

and  (5.9c))  develops  a  logarithmic  singularity,  indicating 
breakdown  of  the  solution  and  the  end  of  the  explosion  stage. 

S.3  The  Outer  Region  OR 

This  region  remains  essentially  stationary,  and  hence  plays  no 
role  during  the  explosion  stage.  For  the  sake  of  completeness,  we 
give  below  the  asymptotic  form  of  the  outer  solution  as  x  ■+  0; 
these  expressions  are  determined  by  combining  (4.31)  and  (4.34): 

(5.32a)  T  ~  1  9  1  [-  <  ( 2 Y- 1  )  / />  Inx  -  «n(A1<x/)]  +  ..., 
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-J.--l.-V-i 


■  *«  -  *  "  •  "  '  .nvn  --wir*  \-v  *v  ■vav 


(5.32b) 

p  -v  1  +  0  1 

(5.32c) 

u  ~  Q  *  Aj<x 

-  £n  ( A 

Similar  expressions  can  be  written  tor  p  and  Y.  The  important 


created  by  the  receding  boundary  layer.  In  this  region  T  must 
vary,  at  leading  order ,  from  the  outer  value  1  to  the  inner  value 
1/ ( 1-a) . 

5.4  The  Intermediate  Region  IR 

This  region,  because  of  its  passive  character ,  will  only  be 
described  very  briefly.  It  is  governed  by  the  variables  a  and  X, 
where  X  is  defined  by 


(5.33)  x  =  e  0X. 


Matching  with  the  neighboring  regions  is  carried  out  at  fixed  0,  by 
setting 


X  =  -9  *  *nx 


as  one  approaches  the  outer  region,  and 
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as  the  boundary  layer  is  approached.  Therefore  the  range  of  X  is 

(5.34)  O  <  X  <  <t(  1-p)  ; 

recall  that  u  Mas  defined  in  (4.23).  From  (5.20)  one  can  easily 
conclude  that  leading-order  matching  with  Lg  requires 

(5.35a)  T  ~  1  /  (1  <7 )  +  ...  ,  as  X  <r(l-p). 


with  analogous  expressions  for  p  and  Y,  while 
(5.35b)  p  -v  1,  u  •*»  e-®XO'-l)/y  £(y-i)/y2>AiaCT_ 


Therefore  the  solution  is  sought  in  the  form 

(5.36a)  $  -v  §°  +■  ...,  for  $  ■=  T,  p,  p  and  Y, 


and 


(5.36b) 


u  'v 


e-ex(y-n  /  y  uo 


+  . . 


In  the  (X,ct)  variables  the  full  equations  (2.1)  read 


P q  ~  (up)  x  e  X)  =  0, 

p[uCT  -  uuy  e-0(cr-X)J  _  (i/y>p^  B-©(cr-X)  _ 

p[to  -  uTx  e“®<cr'"X>]  -  f(y-l)/y>[pa  -  upx  .  w, 

pFy„  -  uyx  e-0«*-x>i  = 


—  ( 1  / fl)  w 
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P  -  PT , 

where  W  retains  the  definition  (5.3).  Clearly,  the  solution  is 
stationary,  i.e.,  independent  of  ct  to  all  algebraic  orders. 
Specifically,  the  reaction  term  W  is  exponentially  small  since  one 
expects  T  <  1/(1-ct).  The  1 eadi ng-order  terms  can  then  be 
determined  simply  by  appealing  to  the  matching  conditions  <5.35), 
and  one  finds  that 

<5. 37a)  T°  =  l/Cl-(l-p)_1X] , 

with  analogous  expressions  for  p  and  Y,  while 
(5.37b)  p°  =  1,  u°  =  i  <  y-1 ) /y2yf\la  (I -y)-1X. 

It  is  a  simple  matter  to  check  that  the  above  solution  also  matches 
with  the  outer  expansions  (5.32)  as  X  ->  0. 

The  analysis  of  the  Type-B  explosion  is  thus  complete. 

6.  The  Type-I  Explosion 

In  this  section  details  are  largely  omitted  and  emphasis  is  on 
the  results,  since  the  treatment  follows  closely  the  Type-B 
analysis  just  concluded. 

6.1  The  Induction  8tege 

Figures  7(a-d)  display  the  numerical  solution  of  the  induction 
problem.  The  graphs  are  self-explanatory.  Similarity  with  Figures 
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2  (a-d)  is  obvious,  but  two  points  of  contrast  are  noteworthy. 

Pifst ,  the  temperature  pro-file  now  has  a  rounded  peak.  Second,  the 
boundary  layer  is  thicker;  this  can  be  seen  more  clearly  in  Figure  0, 
where  T j (x ,t ) /T j (0,t)  is  plotted  at  the  last  successful  time  step 
for  each  of  the  two  cases. 

The  boundary  layer  retains  the  form  (4.8)  and  a  two-sublayer 
structure  emerges  once  again.  The  expansions  are 

(A.  la)  T  ~  1  +  0~1[~!n<)'r)  -  rinr  f(/-l)/y>A1  + 

r  f-AlS2  +  A2>  +  . . . ]  + 

(6.1b)  p  ~  1  +  0-1 [-in (Bj r)  -  rinr  CCZ+D/XlAj  + 

r  c-Axs2  +  c < y+i ) / <y-i > } a2  +  a 1/yi  +  ...]  +  ..., 

(6.1c)  u  ~  [rinr  f-2Ajs/y>  + 

r  C2A2/(/-l)  -  Ai/yis  +  ...]  +  ... 

in  L j ,  and 


(6.2a) 

T 

~  1 

+ 

e  1  [-In  ( /r)  -  Jnd+A^2)  +  ...]  +  ..., 

(6.2b) 

P 

1 

+ 

B-1  [—In  (Bj  r )  -  ind+A^C2)  +  ---]  +  ..., 

<6. 2c) 

u 

-V  9~ 

■1 

[f1/2*nr  f-(2/y>A1C>  +  r1/2  C-(2//)A1C  Ind+Aj?2) 

+  C2A2/<y-l)  -  Aj/XKJ  +  ...]  +  ..., 


in  Lg-.  The  coefficients  A^  ,  A->  and  Bj  are  to  be  determined  by 
matching  with  the  outer  solution  as  before.  The  spatial  coordinate 
C  in  l_£  is  defined  by 

1/2 


(6.3) 


<  =  x/c 


1  /? 

implying  that  the  boundary  layer  is  now  0(r  )  thick,  and  hence 

thicker  than  the  Type-B  layer.  The  -finding  o-f  Figure 

8  is  thus  confirmed. 

It  turns  out  further  that  the  L^-solution  is  uniformly  valid  all 
the  way  to  <  =  0,  so  that  the  interior  sublayer  is,  in  fact, 
superfluous. 

For  smooth  merging  with  the  boundary  layer  the  outer, 
numerically  computed  solution  is  required  to  have  the  asymptotic 
form 


(6.4a) 

T  ~  1 

+  B~1 

C-2£nx 

-  £n(/A1) 

(6.4b) 

p  -v  1 

+  e_1 

C— 2£nx 

-  £n (AjB ^  > 

(6.4c)  u  ~  0  1  C -4  ( A  j  /  y)  x  #nx  +  02  <  Aj  /  Y)  InAj 

+  2A2/O'— 1)  —  Aj/y>x  +  ...3  +  ...,  as  x  ■*  O. 

This  behavior  was  confirmed,  and  the  constants  A^,  A-?  and 
computed,  by  comparing  the  numerical  solution  with  the  above 
expansions.  The  remaining  variables  Pj  and  can  be  determined 
as  before,  by  appealing  to  the  first  member  of  (3.2c),  and 
(3. 2d) . 

6.2  The  Explosion  Stags 

The  analysis  proceeds  as  in  section  5.2.  The  appropriate 
coordinates  are  cr  and  C,  and  the  requirement  of  matching  with  Lj 
is  replaced  by  the  condition  of  regularity  at  (  =  O.  The  solution 


turns  out  to  be 


(6.5a) 


T 


+  .  . 


•  » 
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(6. 5a > 

T  ~  T00  + 

8  1 

[T01-(l-,)-2lna+A1(2)]  +  ..., 

(6.5b) 

P  ~  Too  + 

8~X 

[T01—  *  1— ct) _2fnCl+A1 C2>  +  P01/ ( l~a) ] 

+  .  .  .  , 

(6.5c)  u  ~  B~0a/2  Ccr/  C>'(  1-ct)  >  ]  +  ... 

( 6 . 5d )  p  -v  1  +  0-1  p01  +  .  .., 

(6.5*)  Y  -v  Y00  +  (0a/)~1  [^n(b1/y)-T01 

-Kl-CT)-2£n(:i+A1C2>]  + 

where  the  constant  bj  appearing  in  (6.5e)  was  defined  in  (4.3B). 

The  doubly  subscripted  quantities  correspond  to  the  spatially 
homogeneous  explosion,  and  were  introduced  in  (5.9).  The  explosion 
stage  peaks  just  as  it  did  for  Type-B,  and  the  remarks  at  the  end 
of  section  5.2  remain  valid.  Finally,  the  IR-analysis  of  section 
5.4  carries  over,  with  obvious  modifications. 


7.  Concluding  Remarks 


The  spatial  structure  and  temporal  evolution  of  a  localized 
thermal  explosion  in  a  confined  gas  have  been  described 
mathematically.  Localization,  rather  than  a  spatially  uniform 
explosion,  occurs  as  a  result  of  system  nonhomogeneities,  here 
modelled  by  a  slightly  nonuniform  initial  temperature.  Attention 
is  confined  to  what  may  be  called  the  f ast-reaction  limit, 
character i zed  by  the  initial  induction  time  of  the  reaction  being 
comparable  to  the  initial  acoustic  time  across  the  vessel,  so  that 
diffusion  plays  no  role.  This  limit  can  be  achieved  if  the  initial 
temperature  of  the  unreacted  gas  has  been  raised  to  a  sufficiently 
high  level,  perhaps  by  the  passage  of  a  strong  shock.  By  contrast, 
the  si ow-reacti on  limit  would  correspond  to  the  induction  time  and 
the  conduction  time  being  of  the  same  order.  The  latter  problem 
was  the  subject  of  Poland  and  Kassoy's  investigation  C63. 

The  explosion  is  shown  to  develop  in  two  distinct  stages.  The 
first  stage  is  induction,  characterized  by  small  perturbations 
about  a  spatially  uniform  state,  where  the  primary  interaction  is 
between  linearized  acoustics  and  weak  but  nonlinear  chemical 
heating.  Chemical  amplification  leads  to  localized  thermal 
runaway,  or  blowup  of  the  perturbations,  at  a  time  and  location 
determined  by  the  initial  and  boundary  conditions.  The  spatial 
structure  at  blowup  is  sel f -si mi l ar ,  differing  slightly  depending 
upon  whether  the  runaway  site  is  at  the  boundary  or  in  the  interior 
of  the  domain. 

Induction  is  followed  by  explosion,  characterized  by  0(1) 
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variations  in  the  state  of  the  gas.  The  characteristic  chemical 
time  plunges  dramati cal  1 y .  The  acoustic  time  drops  as  well,  but 
not  nearly  in  the  same  proportion,  so  that  explosion  is  dominated 
by  chemical  heating.  There  is  no  time  -for  expansion,  with  the 
result  that  changes  in  the  velocity  and  density  fields  are 
negligible.  Thus  the  gas  explodes  locally  at  essentially  constant- 
volume  conditions,  with  little  change  in  the  spatial  structure  that 
it  inherited  at  runaway.  (Analysis  in  the  Appendix  shows  that  if 
thermal  expansion  is  admitted,  the  corresponding  spatial  structure 
is  necessarily  singular.)  The  explosion  stage  ends  when  temperature 
and  pressure  within  the  explosion  have  peaked,  the  final  values 
being  exactly  the  same,  to  leading  order ,  as  in  the  spatially 
homogeneous  case.  The  subsequent  expansion  of  the  hot,  highly 
compressed  gas,  and  the  eventual  development  of  a  blast  wave,  are 
currently  under  study. 
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Appsndlx 

The  setup  (4.8),  employed  in  the  text  -for  analyzing  the  spatial 
structure  of  blowup,  is  based  on  the  numerical  observation  that 
both  Tj  and  p^  exhibit  identical,  -inr  behavior  as  r  ■*  0. 

This  observation,  found  to  hold  for  all  the  numerical  runs 
undertaken,  implies  that  blowup  is  a  constant— vol ume  process,  since 
density  perturbation  Pj  =  Pi  ~  Tj  remains  bounded. 

Let  us  now  consider  the  possibility  that  for  some  initial 
conditions,  blowup  lies  partway  between  a  constant-volume  and  a 
constant-pressure  process,  and  ask  whether  a  self -similar  structure 
consistent  with  this  nation  exists.  Accordingly,  we  replace  (4.B) 
by 

(A. la) 

(A. lb) 

(A. lc) 

where 

O  <  A  <  1. 

The  case  A  =»  0  corresponds  to  a  constant-pressure  situation,  and 
A  «■  1  to  the  constant-vol ume  case  already  discussed.  Substitution 
into  (4.7),  followed  by  some  rearrangement,  yields  the  leading- 
order  structure  equations 


Tj  ~  -inr  +  fQ(s)  +  ..., 
Pj  ~  -Afnr  +  gQ(s)  +  . .., 
U  1  -V  f>Q  (s)  +  .  .  .  , 
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<A.2a)  s<l-s2)f0'  +  Cl-s2-(<y-l)/y>AD  -  (l-^lexp  fQ, 

(A. 2b)  (l-s2)gQ’  =  sCA-/exp  f , 

(A. 2c)  y<  l-s2)h0'  =>  /exp  fQ  -  A. 

The  transformation 

(A. 3)  -f  0  =  -fnF 

reduces  (A. 2a)  to  the  linear  equation 

(A. 4)  s  < 1-s2) F '  -  Cl-s2-t(y-l)/y>A3F  =  l-/s2. 

□nee  F  is  known,  fQ,  9q  and  hQ  can  be  computed  sequentially  from 
(A. 3)  and  (A.2b,c). 

Equation  (A. 4)  has  singular  points  at  s  =  O  and  1.  It  can  be 
shown  that  in  general  integration  can  remove  at  most  one 
singularity,  thereby  yielding  solutions  which  are  singular 
either  at  0  or  at  1 .  Such  solutions  can  evolve  only  from  very 
special,  singular  initial  conditions,  and  are  therefore 
unacceptable  if  the  initial  data  are  smooth.  The  only  regular 
solution  is  the  constant 

f  =  y, 

which  requires 

A  =  1, 

correspondi ng  to  the  constant-volume  blowup  already  discussed. 
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Figure  8 

Plots  of  T i (x , t ) /T j <0 , t )  at  the  last  integration  step  for 
(a)  Type  B,  and  (b)  Type  I. 
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Mi  A  list  rad 


Experimental  observations  of  ignition  in  premixed  gaseous 
reactants  indicate  that  perfectly  homogeneous  initiation  is 
practically  unrealizable.  Instead,  combustion  first  sets  in,  as  a 
rule,  at  small,  discrete  sites  where  inherent  inhomogeneities  cause 
chemical  activity  to  proceed  preierentially  and  lead  to  localized 
explosions.  Combustion  waves  propagating  away  from  these  "hot 
spots"  or  "reaction  centers"  eventually  envelope  the  remaining 
bull  . 

This  study  examines  the  spatial  structure  and  temporal  evolution 
of  a  hot  spot  for  a  model  involving  Arrhenius  kinetics.  The  hot 
spot,  character l z ed  by  peaks  in  pressure  and  temperature  with 
little  diminution  in  local  density,  is  shown  to  have  one  of  two 
possible  self-similar  structures.  The  analysis  employs  a 
combination  of  asymptotics  and  numerics,  and  terminates  when 
pressure  and  temperature  in  the  explosion  have  peaked. 
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